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We study BPS vortex configurations in three dimensional U (N) Yang-Mills theories with Chern- 
Simons interaction coupled to scalar fields carrying fiavor. We consider two kind of configurations: 
local vortices (when the number of flavors Nf — N), and semi-local vortices (when Nf > A''). In 
both cases we carefully analyze the electric and magnetic properties and present explicit numerical 

■ solutions. 
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■ I. INTRODUCTION 

D 
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Non-Abelian vortices may find their application in a variety of problems ranging from particle physics and cosmology 
(e.g. confinement, supersymmetric and supergravity models, hot or dense QCD, cosmic strings) to condensed matter 
physics (e.g. quantum Hall effect). Of particular interest are those vortices solving first order BPS equations, which 
correspond to the saturation of the Bogomolny bound for the mass and naturally arise in supersymmetric theories 
' (see refs. l]-([3] for reviews with complete lists of references). 
Qh' Non-Abelian BPS equations for vortices have been analyzed both for Yang-Mills-Higgs and for Chern-Simons- 

^ Higgs models. In the former case the gauge and the (fourth order) Higgs potential coupling constants have 

to be related in order to pass from the second order equations of motion to a first order BPS system. When the 
Chern-Simons term dictates the dynamics of the gauge field, one is forced to choose a sixth order Higg s potential in 
order to find a Bogomolny bound for the vortex mass. Again, coupling constants should be related jl6l|-[l7t. pd| . 
, The origin of these requirements can be also understood in the framework of supersymmetry: they are necessary 
conditions for the existence of an A/" = 2 supersymmetry extension of the bosonic models. In this context, the first- 
order BPS equations arise studying the supersymmetry algebra and looking for supersymmetric states. The resulting 
. selfdual and anti-selfdual solutions break 1/2 of the original supersymmetry [l8l|-[20l|. 

The mixed case of Yang-Mills-Chern-Simons (YMCSS) vortices was also recently discussed [llj-jl^. As in the 
abelian case [l^, in order to have a Bogomolny bound and first order BPS equations the coexistence of the two terms 
' giving dynamics to the gauge field requires a careful choice of the number and type of scalars. Moreover, the Gauss 
, law, through which the Chern-Simons term enters into the energy, is no more an algebraic equation for as in the 
case when the Yang-Mills term is absent, but a second order differential equation that should be taken into account 
together with the first order BPS system. 

It is the purpose of this work to construct explicit BPS vortex solutions for the the YMCSS model thus completing 
5h ■ the analysis presented in [22] where the low-energy vortex dynamics was the main aim of study. We shall take U{N) as 
gauge group and include Nf scalars in the fundamental representation and one real scalar in the adjoint representation 
of U{N). We shall consider both cases: when the number of flavors Nf is equal to the number of colors N and also 
when Nf > N, in which case vortices become semi- local [2^.jl^.[2H-[26j. 
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II. THE MODEL AND THE BPS EQUATIONS 

We shall consider the d = 2 + 1 dimensional U(N) Yang-Mills-Chern-Simons model discussed in ^] with dynamics 
governed by the Lagrangian 

£ = - -i^TrF^.F'^'' - -^Tie^^'^P ( A^d^Ap - ^A^A.Ap) + \tt{V^ 



Here are A'^j scalars with i the flavor index (i = 1, 2, . . . , A^/). Each transforms in the fundamental representation 
of the gauge group U{N) and is a real scalar in the adjoint. We shall first discuss the Nf = N case (local vortices) 
and then extend the analysis to Nf > N (semilocal vortices) Whenever it does not lead to confusion summation over 

flavors is implicit. The gauge fleld A^ takes values in the Lie algebra of U{N), A^^ = A^t^, where t"^ = (1"^)°^^ are 
the U{N) generators {A = 1, . . . , - 1; a, 6 = 1, ... TV) with normahzation Tr t^i^ = 5^^/2. The curvature and 
covariant derivatives are defined as 

Ffj.f — dpAn — dyAp — i [A^i, A^] 
'Dpcj) = dp(j)-i [Ap,4>] 

T^fj^Qi = df^qi - iApQi (2) 

The Chern-Simons coefficient k must be an integer for > 1. 

The masses rrii of the fundamental scalars break the flavor symmetry to U{1)^~^. There exists a fully broken Higgs 
phase in which the scalars take the following expectation values 



<lt..c^St\lv^ + ^ <f>tc = S'^'m, <_ = 0. (3) 

In this vacuum, the U{N) gauge symmetry and the U{1)^^^ flavour symmetry are spontaneously broken. There 
remains only a diagonal symmetry, U{N) x U{l)'jf^^ f^(l)dtag^- This corresponds to the combined action of a 
gauge group element € U{N) and a flavour transformation Vij ~ %e*"^ S U{1)^~'^ (with — 1 independent 
parameters aj): 

Cac - U^^'..cVJ^ - U'^'^lt {U-'f (4) 

where U = . 

Let us note that in the ^ 00 limit where the Yang-Mills term and the kinetic energy term for the 4> fleld can 
be discarded, the adjoint fleld may be eliminated from the Lagrangian, and for nii = one ends up with the sixth 
order potential which allows one to flnd BPS equations for the pure CS-Higgs system both in the Abelian p^-(l7j 
and non-Abelian JJ^ cases, 

hm V[q,cl>,m,=0]=^-^(\q,\'-v'Y\q.,\\ (5) 

The energy associated with Lagrangian ([T]) can be constructed from Tqq, the time-time component of the energy 
momentum tensor T^i,. It takes the form 



E = I (fx Too = / fx 



^-Tr {El + S2) + iTr {{Vo^f + {V^<i>f) + \Doq^\^ + \V^q^\^ 

4. 2 ^ ( \ ^4' ? 

+ 9i (0 - mi) q, + — Tr I q,q\ ~ 1^ - ^ 

where Ea ^ Foa, B ^ F12. 

Following the Bogomolny procedure of square completion a lower energy bound can be obtained [l^l 

E> 2T:nv'^ +^Q,m,. (7) 



^ 2 
,2 ' 



(6) 
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Here n G Z is the topological charge of the configuration, which corresponds to the topological degree associated to 
the qi component that carries the winding, 



n = — Tr / (fxB 

and Qi are the conserved Noether charges associated with the residual U{1)^~''^ flavor symmetry, 

Qi = i J (fx {qlVoqi - (Voqi)^ qi^ ■ 
Using Gauss law one can find the typical Chern-Simons term connection between charge and flux, 

i 

with $ the magnetic flux, 

$ = ^Tr / (fxB = —n 
e J e 

The bound ([7]) is saturated whenever the following BPS first-order equations hold 

Viq^ ± il?29i = 

, i K(j) 
B±-iq,qj-^-v') = 
Z Zn 

T^oqi T mi)qi = 0. 

It should be signaled that in obtaining the bound, it is necessary to use Gauss' law 



K I 



{V^qi) q\ ~ q., [Voqi) 



1 



-T>aEa 



= 0. 



which should then be considered together with eas. p^ - p^ when looking for explicit vortex solutions. 
At the bound, the energy of configurations can be identified with the BPS soliton mass. 



M = 



(8) 

(9) 
(10) 

(11) 

(12) 

(13) 
(14) 
(15) 
(16) 

(17) 



(18) 



In what follows we choose the upper sign in eas. p2|) - ^6l) . which corresponds to a non- negative winding. Of course, 
the opposite choice is equally treatable. 
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III. THE VORTEX ANSATZ 

Starting from the trivial vacuum, a winding can be introduced through a singular gauge transformation generated 



VL{ip) = diag [1, 1, e™"^] = e'^diag[e-*-&'^, e"*^"^, e*" 



(19) 



We have written the formula above, so as to emphasize that combines an C/(l) element with an element of Zat, the 
center of SU{N). Then, a configuration of the form 



^{(p)qvac 



(20) 



with q^ac the trivial vacuum ([3]) will lead to a topologically nontrivial but singular (at the origin) string configuration. 
To avoid the singularity the natural ansatz for a regular vortex should be 



q = diag [m,^, ...,T]Ne'"'^qN{p)] 



(21) 
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with with "ql = + Kmi/2'K and QNip) vanishing at p = 0. Then, (I2ip should be supplemented with consistent ansatz 
for the remaining fields, 

(f> = diag[TOi,TO2, ■•■,mAr + /lAr(p)] (22) 
= diag[0,0,...,n-ajv(p)] (23) 
Ap = 0, (24) 

The complete set of appropriate boundary conditions ensuring finite energy is: 

aN{0) = n, a7v(oo)=0, gjv(O) = (25) 
qjv((X3) = 1, /i^(oo) = 0. (26) 

Concerning Aq, equations p4|) and p5|) require 

[^o>]=0 5p(Ao + 0) = O. (27) 

This suggests 

Ao = -(^ + C with [C, 0] = 0, (28) 

with C determined by (fTB|): 

(</. + Ao)°'g,'' = m,gf ^ C"'^ = ,5''Vfc. (29) 
Unless all masses vanish, one cannot set Aq = —(p. The above equation fixes Aq in our ansatz: 

= diag[0,0,...,-/i7v(p)]. (30) 



IV. THE BPS VORTEX SOLUTION 

Plugging ansatz ([211)- dSJ) into equations ^ and (d?]) gives 



pdpqN-aNqN^O (31) 

1 6^ / /S^ \ 

- dpQN - y ('7^9^ - -^hN - 77w j = (32) 



dpttN - hNvff'lN + ^ (dpiN + ^ dpIlN^ ^ 0, 



(33) 



Note that the Gauss Law constraint is at the origin of the second-order derivative in /ijv in p3|) . It will be convenient 
to define 



„2„2 



and 



a(r) = a7v(p) g(r) (7Ar(p) /i(r) = lhN{p) (35) 
where t = \fPp. One can use ((3T|) to eliminate dpa^ in (133^ so that equations can be recasted in the form: 



T (q^ -h~l) (36) 



da 

— ^ - qa (37) 
dr T 

dh 

— = (38) 

^ = 2/iq2 - a (g2 _ _ 1) _ i (39) 
dr ^ ' T 
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where a = /3j'^. The boundary conditions imposed by finite energy read now 

a(0) = n , a{oo) — 
g(0) = , <7((X3) = 1 
h{oo) = 0. 

Concerning the behavior of h at the origin, it should go to a finite constant. 
Using (fT8|) . the vortex mass for our ansatz takes the form 

M = 2nv'^n + Qnitin — 27r77^n. 

The BPS vortex mass is solely determined by the topological charge, the role of 77^ being just that of a scale. 
Our result is of course consistent with the re-parametrized form of the energy 



(40) 



(41) 



E = 



2n 



rdr 



{■^ (h'f + ^ («')' + ^q'h^ + ^^aV + 2(3 {q'f +(3{q^-h- if^ , (42) 



which after some algebra and integration by parts can be written as follows: 



27r 



1 Tdr (2/3 




2 




L r J 





h-1] 



- {Th')' ~2hq^±a- 

T T 



(43) 



The upper sign corresponds to our non-negative winding ansatz. The first three terms in the integral vanish as they 
are readily identified with the Bogomolny equations, while the last term gives the expected contribution to the vortex 
mass M = 2'iTrij^n. 

Let us note that for niM — 0, eqs. (|36l) - (|39|l reduce to the abelian case discussed in [2^]. It is indeed typically 
observed that the non-abelian Z^r-vortex equations of a model reduce to the abelian equations when the coupling 
constants of the U{1)- and the S'f7(iV)-gauge groups are set equal (this choice has been made implicitly here by 
working with the gauge group U{N)). In the case rriiq ^ 0, the only modification with respect to the mjv = case 

arises through the parameter a = {2tiv'^ -\- ktun) ■ Hence we can obtain the profile functions of any general 
{mjv ^ 0, k} vortex from a {toat = 0, k'} solution by setting k' — k/ -y/l -|- nmpf /2ttv'^ . However, the behavior of the 
physical observables depends upon the value of m^. If the Chern-Simons term is absent {k = 0) the Gauss law is 
satisfied hy h — and our ansatz reduces to the well-honored Abrikosov-Nielsen-Olesen vortex |27|-(28j. 



To obtain numerical solutions to the BPS equations ([36|) -(|39 p we have used a relaxation method, selecting the following 
four boundary conditions 

a(0) = 1, q{oo) = l, h{oo)^u{oo)^0. (44) 

Given our ansatz, the magnetic field B and the electric field £ can be defined as 

B = TrFi2 , £ = Tri^^op- (45) 

They are depicted in Figure [T] where we have set e^/v^ = 1. It is interesting to observe the behaviour of the magnetic 
field in the case mpf ^ 0. As k is increased, the magnitude of the magnetic field at the origin initially increases. At 
large enough k, however, the i3-field starts to decrease at the origin and begins to show a characteristic bump, as 
encountered in the Abelian case [2^. Any further increase in n amplifies the size of the bump. The electric field is 
considerably smaller than the magnetic field for small k, but it also increases as the CS-coupling becomes important. 
As in the m^r = case, the ratio i£'max('«)/'Smax(«^) increases linearly with n for small k, and eventually tends to a 
constant. 

We may explain this behavior by observing that we have q;(k) cx V k when tojv — and a{K) cx k for 
K S> 2iTv'^ /iTLN when ra^ 7^ 0. The limit k 3> 2t:v'^ /mM is where the effects of non-zero ttim become noticeable. The 
profile functions a, h and q thus depend more sensitively on k when the mass is zero. Furthermore, when computing 
the electric and magnetic fields from the functions a and /i, their behavior is not explicitly dependent on k for m^r = 0, 
whereas we have explicit dependence of order 0{k) in B and of order ©(-^/k) in £ for mjsi 7^ in the same limit. In 
the zero-mass case [23|, the magnetic field starts to exhibit the typical doughnut-shape as k, is increased, at the same 
time as its overall strength decreases. The electric field develops the characteristic bump and also becomes weaker. 
As a consequence, the ratio fmax/S,nax, which vanishes as k — > 0, approaches a finite constant for k 3> 2'kv'^ jmj^. 
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Figure 1: The magnetic field B and electric field E of the [n = 1, niN /v^ = 1, k} vortex are shown for k — Q, 10, 60, 100. The 
K = line corresponds to the abelian Nielsen-Olesen vortex, which exhibits no electric field. 



For TO jv 7^ 0, on the other hand, the explicit K-dependence of the electric and magnetic fields counter-acts the trend 
of a and h becoming smaller for larger k. The maxima of the two fields now increasing for a larger range of n and 
do not go to zero as k — > oo. As a result, the ratio fmax(«^)/'Bmax('«) still tends to a finite constant asymptoti- 
cally, but it reaches it more slowly as m^r becomes larger. It is also worth noting that the radius of the bump in 
the electric field linearly increases with k, when tojv = 0, whereas it approaches a constant when tojv = (see Figure[T|). 

Let us end this section by computing the angular momentum for the vortex solution using the formula 

J = / (fxSijX^Toj^ / (fxTo^p (46) 



(which is actually independent of the ansatz). Given Lagrangian ([T]), one has 

Tov = ^pTr {Ep B) + ^Tr (Doc^D^cj)) + ((i^o^a)^ (D^Qa) + {D^Qa)^ (i?o<Za)) 
2 

= —h'^a'^ + 2ri%hN Qn (47) 
Using Gauss' law, which for our ansatz reads 

V% h^ql^^ (^h% + i h'^^ + -^{a'^/p) (48) 
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we can write 



So, we have for J 



To^^——{paMh'^) + - (49) 

p dp An p dp 



2 K 
J = 27r ( -i^paN h'j^ + —a% 
4 7r 



(50) 



or, finally 



J^-^n^ (51) 

Being k an integer, also the angular momentum is quantized at the classical level. Note that in view of eq. (|10p which 
in the present case reduces to Qn = nn we see that the angular momentum can be written in terms of the square of 
the charge (in contrast with the pure Yang-Mills case in which it is proportional to the charge [oj-fiojl. 

V. SEMI-LOCAL VORTICES 

Unlike Abrikosov-Nielsen-Olesen vortices and their non-abelian extensions, the radius of semilocal vortex is not 
fixed but it becomes a parameter. This kind of vortices emerge when Nf > N, that is, when there are = Nf — N 
additional fundamental scalars {pe — qn+c}, e = 1, ...,-/V/ — iV in comparison with the local vortices arising in the 
N — Nf case. 

We then start from Lagrangian ([1]) now with Nf > N and consider minimization of the potential in the general 
case nii ^ rrij ^ i j, with i,j < N. The mass term in the Nf > N Lagrangian 

^ qj'^ {r'^ - S'^'m^fql (52) 

can only be made to vanish for N non-zero fundamental scalars, since we may only pick N of the diagonal entries in 
(j> to cancel the S'^^nii terms. Then, in order to minimize the potential the remaining qi need to vanish. Without loss 
of generality, we may choose the {pe — qN+e\ to vanish. These fields then lie in the unbroken vacuum. Spontaneous 
symmetry breaking of U{N)g x U{1)^^ ^ occurs as before for the original TV-sector while the additional fields exhibit 
invariance under the more general transformation 

p%,, exp(*ae)C/"'p^vac, (53) 

with unconstrained {ag}. These fields must be topologically trivial. If we adopt the previous ansatz for the original 
fields, equation which as we explain below is still valid for Nf > N requires all components of to vanish 

identically except for a = N. This suggests the following ansatz for the Pe 

p- = riNS^^Ur). (54) 

Furthermore, it is required by that the masses of any additional (non-trivial) scalar pe be equal to the mass of 
the field that carries the winding, 

TOg = mN if Ceir) ^ 0. (55) 

Equation ^ for the energy in the Nf — N case is still valid for Nf > N and so is the bound ([7]) and the 
BPS equations P^ -(|16 p . Concerning the axially symmetric ansatz, it consist of the one proposed in the local case, 
eas. (PT|) - (l^^ . augmented with (|54p for the extra scalars. Inserting the ansatz in the BPS equations one now obtains 

^=r{q' + \a'-h-l) (56) 

— ^-qa 57 
dr T 

^ = i(a-n)ee, e^l,...,Nf-N, (58) 
dr T 

dh 



^ = 2k(,» + 161") - a(,= + If. |» - - 1) - 1 (60) 
dr T 
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with l^el = Eefke- Equations ^ and §^ can be used to solve for the profile functions {^e} [2j 



Ce(T) = X. 



(61) 



with Xe € C arbitrary complex constants that parametrize the solutions. Of course, if we set all the Xe parameters 
to zero the system (|56l) - ([5n)) coincides with (P5 )) - (|39p and the semi-local vortices become ordinary local ones. 

The energy of the semi-local vortex is 

Es=^^J rdr (^A {h'f + A ^af + ±h^ [q^ + + 2^aV + 2^ia - nf 

+2/3 ({q'f + \i'f) +P(q^ + -h-lfy (62) 

As in the local case, it is straightforward to show that this can be written as: 



27r 
72" 



j rdr (2/3 


W T ^1 


2 

+ 2/3 


^, (a - n)Ce 

Se + 


2 


'a' 




L T J 




T 




T 



+ 2/3r-ia' 



(63) 



the upper sign corresponding to the non- negative winding vortex. Energy (j63p reduces to the lower bound in when 
the Bogomolny equations are satisfied, 



27rnu^ + ^ 



(64) 



In the case of semilocal vortices, one can define a parameter x, the complexified size of the vortex through the 
formula 



(65) 



We thus see that, as expected, the vortex mass is x-independent but the behavior of the fields at infinity drastically 
changes with respect to the local case: the fields have a long range power falloff instead of an exponential one. This 
can be seen in Figure [H which shows the numerical solutions to (ISSjl - dBU)) . Indeed, at large distances (p ^ l/e^fp) 
and for very large transverse size of the vortex (x ^ 6%/^) can see how the asymptotic behavior is no more that 
of an exponential falloff but a power one. The analytical asymptotic behavior of the fields aat and q^q are 



.2|n| 



g^r ~ i — -- 



2o2|r 



whereas has the same exponential falloff behavior as in the local vortex case: 



'llN P 



In 



(66) 



(67) 



Let us finally note that the presence of the radius x also reduces the Chern-Simons characteristic bump at the 
origin. 

Concerning the angular momentum for the semi-local vortices, the component Tq^ of the relevant energy-momentum 
tensor component gets just an extra term ATq^ with respect to the local case. 



,11 N 



= + 2— hn{n ~ ajv) ?7n l^e 



(68) 



Figure 2: The magnetic and electric field of an {n = 1, uin /v^ = 1,k = 100} vortex with |x|^ = 0, 2, 5, 50, 100. 



and then 



Using Gauss' Law 



2 

To^ = ±—h'j^ a'j^ ±2if^hN {ql^ + \ie\^) an T '^nri% hn \£,e\^ (69) 



we can write 



P"' J-p^N{p) p) = V% hn {q% + ICel") T ^P"' (70) 

To^ = ±^~{h'N aN p) + ^-^ia%) T^nv^hN (71) 

where the two first term coincide with the local vortex value for Tq,^. Taking into account boundary conditions one 
now has 



J = 27r j To^pdp= -^n^ ±2nr]% j hN\ie\' pdp (72) 



Thus, unlike the case of local vortex, the angular momentum is not simply quantized in terms of K'n? but it depends 
explicitly on the size of the semi-local vortex. 
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VI. DISCUSSION 



A rich spectrum of non-Abelian vortices in theories where the dynamics of the gauge field is governed both by a 
Yang-Mills and a Chern-Simons action was shown to exist in refs.p]|-[23|. where the low-energy vortex dynamics was 
described in terms of a gauged sigma model on the vortex worldline. Although the BPS equations were obtained, 
explicit solutions were not presented and this was precisely the main objective of our work. To this end, we proposed 
an axially symmetric ansatz leading to BPS vortex solutions for a Yang-Mills-Chern-Simons U{N) gauge theory 
coupled to scalars when the number of flavors Nf > N , analyzing the electric and magnetic properties of the local 
{Nf = N) and semi-local {Nf > N) vortices. 

A first interesting feature of the local vortex solutions concerns the localization of the magnetic and electric fields. 
As expected, as the Chern-Simons coefficient k grows, the maximum of magnetic moves away from the vortex center 
and the electric field, also with an annulus shape, starts to develop. Semi-local vortices exhibit a similar behavior 
except that B and £ have a long range power falloff instead of an exponential one (only the real scalar field keeps its 
exponential falloff behavior) . Another difference between local and semilocal vortices concerns the angular momentum 
which is a purely topological object in the former case while it depends on the vortex size in the later semilocal case 

Our investigation started from the Lagrangian proposed in (2^ with the scalar potential and constants chosen so 
as to guarantee the possibility of an A/" = 2 supersymmetric extension, this ensuring the existence of BPS equations 
[l|-[l]. Actually, it would be of interest to investigate the properties of the supersymmetric model and to construct the 
low-energy effective action describing moduli dynamics and analyze its properties both at the classical and quantum 
level, following the approach presented in ref.[l3| for the pure Chern-Simons case. We hope to discuss this issues in 
a future work. 
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